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I. INTRODUCTION 

Heterogeneous magnetic superconducting systems (HMSS) represent a new class of nanos- 
tructures. They are made of ferromagnetic (FM) and superconducting (SC) pieces separated 
by thin layers of insulating oxides. In contrast to the case of a homogeneous ferromagnetic 
superconductor studied during the last two decades, the two order parameters, the magnetiza- 
tion and the SC electron density do not suppress each other llLl^. In HMSS, the interaction 
between the two order parameters is due to the magnetic field created by the magnetic and SC 
textures. Strong interaction of the FM and SC systems not only gives rise to a new class of 
novel phenomena and physical effects, but also shows the important technological promise of 
devices whose transport properties can be easily tuned by comparatively weak magnetic fields. 

The interplay between ferromagnetism and superconductivity has long been the focus of 
studies both experimental and theoretical ||3l|4j,|5j,|6|]. In 1957, Ginzburg pointed out M\ that the 
two phenomena can occur in thin films and wires, due to the small FM induction and relatively 
large SC critical fields. Afterwards, Anderson and Suhl [8] described the conditions for both 
phenomena to appear simultaneously in the bulk. The domain-like magnetic structure in the co- 
existence phase was noted independently by two different groups [9, 10]. Later, Suhl developed 
the detailed version of the Landau-Ginzburg theory of ferromagnetic superconductors [ 11], e.g. 
HoMosSs and ErRli4B4. The first experiment in this context was carried out by Mathias et al. 



1 1211 . The coexistence of ferromagnetism and superconductivity in the bulk has been observed 
recently in the cuprates RuS^GdC^Os-y H13TI and RuS^Gdi+.^Cei-xQ^Oio 1 14, 15] below 
their SC transition temperatures Ti, = 15 — 40 K and T s = 37 K ( for x = 0.2 ) respectively. 
It is possible to avoid the mutual suppression of the FM and SC order parameters by separat- 
ing them in space [1]. Such structures can be made with modern nanofabrication techniques. 
The proximity effect and spin diffusion which suppress both order parameters can be easily 
avoided by growing a thin insulating oxide layer between the FM and SC components. Several 
theoretical studies have proposed different possibilities for realization of HMSS: arrays of mag- 
netic dots on the top of a SC film y],LLa], single magnetic dot of various geometries on top of 
a SC film |l l7 U l8 Ul9 1. ferromagnetic/superconducting bilayers (FSB) and multilayers (FSM) 
[20, 21, 22L|23J, magnetic nanorods embedded into a superconductor [24, 25] and semi-infinite 
HMSS Q. 

In HMSS, the strong interaction between FM and SC components stems from the magnetic 
fields generated by the inhomogeneous magnetization and the supercurrents as well as SC vor- 
tices. SC vortices were widely studied, both experimentally and theoretically, in the literature 

1 



1 2811 . Vortices in SC films were first studied by J. Pearl 12911 . He realized that their current 
decays in space more slowly than those in the bulk. 

Similar models of HMSS were studied theoretically by other groups in different configura- 
tions. Marmorkos et al. investigated the giant vortex state created by a magnetic dot, of size on 
the order of the coherence length L embedded in a SC film by solving the nonlinear Landau- 
Ginzburg equations W17I1 . Kayali IU8LI19I1 and Peeters et al. [30] studied vortex creation by both 
in-plane and out-of-plane magnetized ferromagnets of different shapes on a SC thin film. Sasik 
and Santos et al. considered an array of FM dots on a SC film [31, 32]. They treated the dots 
as magnetic dipoles, ignoring their real geometry and showed that they excite and pin SC vor- 
tices. Carneiro studied interaction between vortices in SC films and magnetic dipole arrays, and 
showed that the arrays with antiferromagnetic order is more effective in vortex creation than the 
arrays with ferromagnetic order [33]. Symmetry violation in SC thin films with regular arrays 
of magnetic dots are reported by two different studies lllo. 13411 . Bulaevskii et al. discussed the 
pinning of vortices inferromagnet- superconductor multilayers [35]. The same group also exam- 
ined the effect of the screening magnetic field on a thick magnetic layer which is placed on top 
of a bulk superconductor [36]. They found that the magnetic domains shrink in the presence of 
the SC film. In the most of theoretical studies, hard magnets are considered. Recently, Helseth 
found that identical vortices can attract each other in the presence of soft magnets 113711 . 



To date, only sub-micron magnetic dots covered by thin SC films have been prepared and 
studied 113 8l l39l l40l l4lL l42l I43L 1441 . The experimental samples of FM-SC hybrid systems were 



prepared by means of electron beam lithography and lift-off techniques 114311 . Both in-plane and 
out-of-plane magnetization was experimentally studied. The dots with magnetization parallel 
to the plane were fabricated from Co, Ni, Fe, Gd-Co and Sm-Co alloys. For the dots with 
magnetization perpendicular to the plane, Co/Pt multilayers were used. The FM dots were de- 
posited on thin SC films made of either Nb or Pb, whose transition temperatures are around 
7 — 10 K. In these experiments, the effect of commensurability on the transport properties ( 
e.g. magnetoresistance oscillations and matching anomalies ) was observed. However, this ef- 
fect is not limited to the magnets interacting with superconductors and was first found many 
years ago by Martinoli and his group [46]. They studied the transport properties of SC films 
with periodically modulated thickness in an external magnetic field. They found oscillations of 
critical current versus magnetic field. Recently, several experimental groups observed commen- 
surability effects caused by a periodic array of magnetic dots or holes Il47l.l48l.l49l.l5 Oil . These 
results confirm that the FM dots create and pin vortices. However, much more interesting and 
promising would be effects specific for the HMSS, which are associated with the violation of 
the time reversal symmetry. Some of them include spontaneous currents in the ground state 
l!24ll5lLl52ll53h . So far only one such effect was experimentally observed: an asymmetry of the 
SC hysteresis in the presence of magnetic dots Il39h . 

In HMSS, the magnetic field induced by inhomogeneous magnetization penetrates into a su- 
perconductor through SC vortices, while the magnetic field generated by the supercurrents and 
SC vortices acts on the magnetic system. The mutual interaction between FM and SC subsys- 
tems offer interesting physical effects such as spontaneous symmetry violation II 1611 . shrinkage 



SC transition temperature shift 12311 



of FM domains and a magnetic domain wall 12 lL l36l 
and Bean-Livingstone-like energy barrier for vortices 

We recently developed a method based on London-Maxwell equations to study theoretical 
realizations of HMSS W27I1 . In this article, we first review this method and its extention to pe- 
riodic systems and finite systems. Next, some of our recent results on a SC film with a FM 



dot grown upon it and domain structures in FSB and further studies in these systems are briefly 
discussed. This review is prepared as a progress report and aims to give reader the basic as- 
pects of HMSS. This article is organized as follows: In the next section, a method to calculate 
inhomogeneous magnetization and supercurrents including SC vortices in the London approx- 
imation is presented with great details. In section II, a FM dot on top of a SC film and vortex 
states in the ground state are discussed. Section III is devoted to the recent theoretical results 
on FSB. We conclude with discussions and summary in section IV. To help reader follow easily 
the theoretical analysis in this chapter, we give details of calculations with Bessel functions and 
series in the appendices. 

H. THEORY 

In both theoretically proposed and experimentally realized HMSS, the magnetic texture in- 
teracts with the superconducting (SC) current. An inhomogeneous magnetization generates 
a magnetic field outside the magnets that in turn generates screening currents in the super- 
conductor, which subsequently change the magnetic field. The problem must be solved self- 
consistently. In the literature, HMSS have been studied through Landau-Ginzburg equations 
11711 and linear London method . Though the SC vortices are treated more accurately in the 
former in which SC electron density n s changes in the vicinity of the vortex core, solving non- 
linear Landau-Ginzburg equations is numerically difficult. In London approach, n s is assumed 
to be constant, and vortices are treated as points. However, London's approximation is sufficient 
when the sizes of all the structures in the problem greatly exceed the coherence length £, and 
offers more analytical insight. For this reason, we have been studying the several realizations 
of HMSS through a method based on London-Maxwell equations. Here we present the method 
to calculate inhomogeneous magnetization and supercurrents including the SC vortices in the 
London approximation. In the next section a method for the most general 3-dimensional HMSS 
is given. In section B, this method is applied to the case of very thin ferromagnetic (FM) and SC 
films. In section C, the method is extended to periodic heterogeneous magnetic superconducting 
systems (PHMSS). In the last section, we show how this method can be applied to semi-infinite 
and finite systems. 



A. Three Dimensional Systems 

The total energy of a stationary FM-SC system reads 15 511 

where B is the magnetic induction, M is the magnetization, n s is the density of SC electrons, m s 
is their effective mass and \ s is their velocity. We assume the SC density n s and the mag- 
netization M to be separated in space. We also assume that the magnetic field B and its 
vector-potential A asymptotically approaches zero at infinity. After the static Maxwell equa- 
tion V x B = — j, and B = V x A are employed, the magnetic field energy can be transformed 
as follows: 



to dV = J-^ dV - <2) 

Although the vector potential enters explicitly in the last equation, it is gauge invariant due to 
current conservation divj = 0. When integrating by parts, we neglect the surface term. This 
approximation is correct if the field, vector potential and current decrease sufficiently fast at 
infinity. The current j can be represented as a sum: j = j 5 + j m of the SC and magnetic currents, 
respectively 1 5611 : 

nfie 2% 

J. = ^— (V<p-— A), (3) 

2m s <|> 

j m = cVx M. (4) 

where cbo = hc/2e is the SC flux quantum. We separately consider the contributions from mag- 
netic and SC currents to the integral ©, starting with the integral: 

^ J ] m -AdV = l - J "(V xM) -AdV. (5) 

Integrating by parts and neglecting the surface term again, we arrive at 

^-Ji m -AdV = ^JM-BdV. (6) 

We have omitted the integral over a distant surface: 

i(nxM)-AdS. (7) 

Such an omission is justified if the magnetization is confined to a limited volume. But for 
infinite magnetic systems it may be wrong even in the simplest problems. 

We next consider the contribution of the SC current j s to the integral ©. In the gauge- 
invariant Eq.©, cp is the phase of the SC carriers wave-function. Note that the phase gradient 
Vcp can be incorporated in A as a gauge transformation. The exception is vortex lines, where cp 
is singular. We use the equation © to express the vector potential A in terms of the supercurrent 
and the phase gradient: 

2tc n s e L 

Plugging Eq.® into Eq.©, we find 



j, • AdV = — j Vcp • j^V - ^3 / f s dV. (9) 
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Since the superconducting current is 

j s = en s v s . (10) 

The last term in Eq.© equals the negative of the kinetic energy and thus exactly compensates 
the kinetic energy in the initial expression for the energy ([l]). Collecting all the remaining terms, 
we obtain the following expression for the total energy: 

f r nfi 2 ,„ N? nfie „ B-M n 



This expression is complete except for a possible surface term for infinite magnetic systems. 
Note that integration volume includes both superconductors and magnets. Eq. dTTT) allows 
us to separate the energy of the vortices, the energy of magnetization and the energy of their 
interaction. Indeed, as we noted earlier, the phase gradient can be ascribed to the contribution 
of vortex lines alone. It can be represented as a sum of independent integrals over distinct 
vortex lines. The vector-potential and the magnetic field can also be presented as a sum of 
magnetization induced and vortex induced parts: A = A m + A v , B = B m + B V , where A^, B^ 
(the index k is either m or v) are determined as solutions of the London-Maxwell equations 
generated by the magnetization and the vortices. The effect of the SC screening of the magnetic 
field due to the magnetization is already included in the vector fields A m and B, n . If such 
separation of fields is applied, then the total energy, dTTT) becomes a sum of terms containing 
vortex contributions alone, magnetic contributions alone and the interaction terms. The purely 
magnetic component can be represented as a non-local quadratic form of the magnetization. The 
purely superconducting part becomes a non-local double integral over the vortex lines. Finally, 
the interaction term may be presented as a double integral over the vortex lines and the volume 
occupied by the magnetization, and is bi-linear in magnetization and vorticity. 

B. Two Dimensional Textures and Vortices 

Below we show a detailed analysis in the case of parallel FM and SC films, both very thin and 
positioned close to each other. Neglecting their thickness, we assume both films to be located 
approximately at z = 0. In some cases we need a higher degree of accuracy. We then introduce 
a small distance d between the films, which in the end approaches zero. Although the thickness 

of each film is assumed to be small, the 2-dimensional densities of super-carriers tv s ' = n s d s 
and magnetization m = Md m remain finite. Here d s is the thickness of the SC film and d m is the 

(2) 

thickness of the FM film. The 3d super-carrier density in the SC film is n s (R) = d(z)nl (r) and 
the 3d magnetization in the FM film is M(R) = 8(z — J)m(r), where r is the two-dimensional 
radius-vector and the z-direction is chosen to be perpendicular to the films. In what follows the 

2d SC density n s is assumed to be a constant and the index (2) is omitted. The energy dTTT) for 
this special case takes the following form: 

f r nfi ,„ n9 nfie „ bm n 9 

where a = A(r,z = 0) and b = B(r,z = 0). The vector potential satisfies the Maxwell-London 
equation, which is derived from the static Maxwell equation V x B = — j, where j is the total 
current density on the surface of the superconductor, and is given by j = (j s + j m )S(z). The 
supercurrent and the magnetic current densities are given in Eqs.QEI). Using B = V x A, the 
Maxwell-London equation reads 

1 271^2/2 € 

V x (V x A) = -~A8(z) + -V<p8(z) +4tiV x (m8(z)). (13) 

k m s c 

2 l 

Here X = *k\jd s is the effective screening length for the SC film, and Xi = ( 4 m ' tC 1 ) z is the 



London penetration depth 15711 . 



According to our general arguments, the term proportional to Vcp in Eq. (IT3t describes 
vortices. A plane vortex characterized by its vorticity n and by the position ro of its center on 
the plane, contributes a singular term to Vcp: 

Vcp (r,r )= /f (r ~" o) , (14) 

|r — ror 

and generates a Pearl vortex vector potential(see Appendix A for details): 

rccbozx(r-ro) rJi(q\r-T \)e-^\ 

A v o(r-r ,z) = — -— , — / — — t dq, (15) 

2n \r — tq\ Jo l+2Ag 

where J\ (x) is Bessel function of the first order. Different vortices contribute independently in 
the vector potential and magnetic field. In the limit of zero film thickness, the usual Coulomb 
gauge, divA = 0, leads to strong singularity in the vector potential, due to the surface currents 
in the SC and FM films, which leads to the discontinuity at z = 0. Therefore, we choose to 
employ another gauge A z = 0. The calculations become simple in the Fourier-representation. 
Following the previous section, we write the Fourier transform of the vector potential Ak as a 
sum Ak = A m k + A v k of independent contributions from magnetization and vortices. Using the 
following definitions of the Fourier transform: 

A k = [ A(r,z)e-'* r - lk ? z d 3 r, (16) 

a q = /a(r,z = 0)e- ! ' qr J 2 r. (17) 

The equation for the magnetic part of the vector-potential reads 

k(k ■ A rak ) - k 2 A mk = ^ - 471* x m q ^ rf , (18) 

where q is the projection of the wave vector k onto the plane of the films: k = k z z + q. An 
arbitrary vector field Vk in wave- vector space can be fixed by its coordinates in a local frame of 
reference formed by the vectors z,q,zx q: 

Vk = V£z + v}q + V i t(zxq). (19) 

The solution to equation (fTSl with A z = is found by taking the inner product of equation ( fT9"l) 
with q, z and zx q, respectively as below: 

I, _ 4nim^ M _aL, 

Amk ~ k z e %v (2U) 



Anim- 



^k—^P 1 ^ (2D 



l 4ni (k z m\ - qm qz 



r ik,d 

W^ ' k 2 



Kk = -—< + ^^ '-<!*■ ■ (22) 



Integration of the latter equation over k z gives the perpendicular component of a q ( m >: 

i 4nXq(m\ + im qz ) ^_ qd 

It follows from Eqs .(127)1 OTT) that a mq = 0. Note that Eq. (ETT) for the parallel component of 

the vector potential AJ nk does not contain any information on the SC film. This component 
corresponds to zero magnetic field outside the FM film. Therefore, it is not essential for our 
problem. The vortex part of the vector potential A v k also has not a z-component since the 



supercurrents flow in the plane. The vortex-induced vector potential is 115 711 



_ 2ity (qxz)F(q) 
Avk " V{l+2Xq) ' ( ] 

where F(q) = Y,j n j e ' qrj is the vortex form-factor; the index j labels the vortices, rij denotes the 
vorticity of the j'th vortex and Tj are coordinates of the vortex centers. The Fourier- transform 
of the vortex-induced vector potential at the surface of the SC film a vq reads 



ityojq x z)F(q) 
q(l +2Xq) 



a vq = 7TT -*r\ ■ ( 25 ) 



We express the energy (TT2l in terms of the fields and vector-potential Fourier-transforms sepa- 
rating the purely magnetic, purely vortex and the interaction parts: 

U = U vv + U mm + U mv . (26) 

The vortex energy U w is the same as it would be in the absence of the FM film: 

nfi 2 f 2% d 2 q 



U v 



f 2% da 

/vcp_ q .(Vcp q --a vq )-^ (27) 



SmJ Y 4 " ^ <|)o ^(2%) 
However, the magnetic energy U„ 



1 f u d2q 

-jm_ q .t> mq — 2 



Umm — ^ / Ml— q'0 mq ^_ 9 (^o) 



contains the screened magnetic field b and therefore differs from its value in the absence of the 
SC film, but it does not depend on the vortex positions. The interaction energy reads 



U mv = — / (V(p)_ q -a raq 

<\m s c J 



2, 



'(271) 

1 f u d2 ° 

2j m -«- K «W) 



m-qb vq — ^. (29) 



Note that only the form-factor F(q) conveys any information about the vortex arrangement. 



C. Periodic Systems 

A periodic heterogeneous magnetic superconducting system (PHMSS) such as a magnetic 
dot array or a periodic domains in ferromagnet- superconductor bilayers can be studied with the 
method described in the previous section. However, it is necessary to modify the equations given 
above for the periodic structures of interest. In this section, we extend the above method to study 
PHMSSs. In doing so, we still assume that PHMSSs are made of very thin magnetic textures 
with the magnetization perpendicular to the plane and SC films. Their energy is calculated over 
the surface of the SC film. We start with the energy of the 2d systems (fT2l . In the plane of 
the PHMSS, the magnetic field b, the magnetization m, the phase gradient Vcp and the vector 
potential a are 2d periodic functions. Therefore, we need to express them in terms of Fourier 
series. For any function /(r), the Fourier expansion is given by 

f(r) = £f G e' G r f G = | / f(r)e-' G - r d 2 r. (30) 

The G s are the reciprocal vectors of the periodic structure of interest, and Si is the elementary 
cell area. We first express a, b, m and Vcp in terms of the Fourier series as in d30l) . then 
substitute them back in (fT2l . Performing the integral over the infinite area of the 2d system 
and using JV( G + G > r d 2 r = &d G qi, we obtain the energy per unit cell u for the 2d PHMSS 
expressed in terms of Fourier components as follows: 

v-, r n^ .,_ . .9 nfie ,,_, . b G m Gn ,„^ 

-=»^i( v »)oi 2 -^ (v,)o -- -V s >- (31) 

The Fourier coefficients of the vector potentials, a m and a v for both the magnetic part and the 
vortices in terms of continuous Fourier vectors are already given in d23l and (l25l) . They can be 
rewritten in terms of reciprocal vectors as 

4nkiG x zm zG 
* mG = 1+2XG ' (32) 

_ ityo(Gxz)F G 
* vG ~ AG(1+2XGY (33) 

Uusing the Fourier coefficients of the magnetic field and the phase gradient: 

b G = z'G x a G , (34) 

(Vcp) G = 2m(G x zF G )/(xG 2 ), (35) 

and replacing the vector potentials in (I3T1) by d32l and d33l . the energy of the PHMSS per unit 
cell is found term by term as 

4o v^ I^gI 2 



4izR 2 ^G(l+2lGy 
m zG F G 



7^1 



( ; x+2XG 



(36) 
(37) 



-2%lT- 
hf l 



G2 l m = G|2 (38) 



c, -+2XG' 



D. Finite Systems 

In the most of theoretical studies, SC subsystem is considered to be infinite size for the sake 
of computational simplicity. Although it is relatively hard to handle SC system's boundaries 
in both Landau-Ginzburg and London equations, several groups have studied finite or semi- 
finite realizations of HMSS. We recently considered semi-infinite HMSS elsewhere [26], by 
benefiting from the ideas of Kogan's study of Pearl vortex near the edge of SC film[58]. Here, 
we propose an alternative method to treat semi-infinite and finite realizations of HMSS. To 
this end, we modify London-Maxwell equation (see Eq. (fT3t ). The boundary condition can 
be incorporated into Maxwell-London equation by using step function, namely supercurrent is 
zero outside the SC system's boundary. Together with step function, Maxwell-London equation 
reads 

V x V x A = — (j,e(r' -r) +j w ), (39) 

c 

where j 5 = ^-rj-Vcp — ^a) and 0(r' — r) is step function that equals or 1 when r is greater 
than the boundary's position r' or otherwise. For semi-infinite system in which 2-d SC film lies 
on x-y plane and its edge is located at x=0, step function is 0(jc). For a finite circular SC disk 
of radius R, step function becomes Q(R — r). Eq. d39l can be solved by similar techniques used 
in the previous sections, however its solution gives rather complicated integral equations. We 
leave the details of the solution to further works. 



III. FM DOT ON SC FILM 

In this section, we review the studies related to the ground state of a superconducting (SC) 
film in the presence of a circular FM dot grown upon it (see (see Fig. [l])). In this case, the mag- 
netization is assumed to be fixed and homogeneous within the dot, and directed perpendicular 
to the SC film. This problem is previously discussed elsewhere 12711 . in which we predicted 
the geometrical pattern formed by vortices in the ground state. Here, we give further details of 
analysis and some new results. The problems we discuss are: i) under what conditions do vor- 
tices appear in the ground state; ii) where they appear, and iii) the magnetic fields and currents 
in these states. As in the previous section, we assume the SC film to be a very thin plane, and 
infinite in the lateral directions. Since the magnetization is confined within the finite dot, no 
integrals over infinitely remote surfaces or contours arise. In the next section, we treat the first 
case using the method described in the previous chapter. 

Let both SC and FM films be infinitely thin and placed at heights z = and z = d, respec- 
tively. The SC film is infinite in lateral directions, whereas the FM film is finite and has the 
shape of circle with radius R (magnetic dot). The 2d magnetization of the magnetic dot is 
m(r) = mzo(R — r)d(z — d), where o(x) is the step function, equal to +1 for x > and for 
x < 0. We first find the vector potential and magnetic field induced by the dot in the presence 
of the SC film, using Eqs. (122123b . The Fourier-component ofd at magnetization necessary for 
this calculation is 

m k = z—A( q Ry k z d , (40) 

q 

where J\(x) is the Bessel function. From d25t with (l40b . the calculations employ Fourier- 
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FIG. 1 : Magnetic dot on a superconducting film. 



transform of the vector potential at the superconductor surface: 

' / X 



iS% 2 lmR 



*mq 



l+2Xq 



(41) 



In the last equation we have replaced e qd by 1 . The Fourier-transform of the vector potential 
, iSn 2 mRJi(qR) \ _ ad 2qk 



mk 



k 2 



+ (e^ d - e - qd ) 



(42) 



l+2Xq 

Though the difference in the round brackets in equation (|42l appears to be small (recall that 
d must be set to zero in the final answer), we cannot neglect it since it implies a finite, not small 
discontinuity in the parallel component of magnetic field at the two film faces. From equation 
(142)1 we can immediately find the Fourier-transforms of the magnetic field components via 



^mk ~~ ?< ?^rak' ^mk 



"*%^mk- 



(43) 



The inverse Fourier-transformation of Eqs. (143142ft gives the magnetic field in real space: 



(44) 



/»oo 

B r m (r,z) = -2%mR / h{qR)Ji{qr)e-^ 

Jo 



2qk 
\+2qk 



sign(z) + sign(z - d) - sign(z) 



qdq, 



(45) 



10 



where sign(z) is the function equal to the sign of its argument. Note that B' m has discontinuities 
at z = and z = d due to surface currents in the SC and FM films respectively; whereas, the 
normal component B z m is continuous. 

Symmetry arguments imply that a vortex, if it appears, must be located at the center of the 
dot. Indeed, for R ^> X, an analytical calculation shows that the central position of the vortex 
provides minimal energy. We have checked numerically that the central position is always 
energy favorable for one vortex. This fact is not trivial since the magnetic field of the dot is 
stronger near its boundary and a violation of symmetry could be naively expected. However, 
the gain of energy due to interaction of the magnetic field generated by the vortex with the 
magnetization of the dot decreases when the vortex approaches the boundary. 

Another interesting problem is the sign of the perpendicular component of the magnetic 
field. The vector potential generated by a vortex is given by Eq. dZ3b with F(q) = 1. The 
perpendicular component of the magnetic field generated by the vortex is 



55, = ^ 



4>0 f°°Jo(qr)e 



-q\z\ 



2n Jo 1 + iXq 



-qdq. 



(46) 



Numerical calculation based on Eqs. (ERIEfoT ) shows that, in the presence of the vortex centered 
at r = 0, B z on the SC film (z = 0) changes sign for some r > R (see Fig. El), but it is negative 
for all r > R without the vortex. The physical explanation of this fact is as follows. The dot 
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FIG. 2: Magnetic field of a dot with and without vortex for R/X = 5 and §o/&ii 2 mR = 0.05. 



itself is an ensemble of parallel magnetic dipoles. Each dipole generates a magnetic field whose 
z-component on the plane passing through the dot has sign opposite to the dipolar moment. 
However, the field exactly over and under the dipole has the same sign as the dipole and is 
strongly singular. The fields from different dipoles compete at r < R, but they have the same 



11 



sign at r > R. The SC current tends to screen the magnetic field of the magnetic dot and have 
the opposite sign. The field generated by a vortex at large distances decays more slowly than 
the screened dipolar field (1/r 3 vs. 1/r 5 ). Thus, the sign of B z is opposite to the magnetization 
at small values of r (but larger than R) and positive at large r. Measurement of the magnetic 
field near the film may serve as a diagnostic tool to detect a SC vortex bound by the dot. To our 
knowledge, so far there are no experimental measurements of this effect. 

The energy of the system in the presence of many vortices with arbitrary vorticities ni can 
be calculated using Eqs. d26H29l ). The appearance of N vortices with arbitrary positions r, in the 
system changes the energy by an amount: 

N j N N 

A N = Y, n ? £ v + 9 L n i n J e w (ry) + Y, ni£mv ( r *) • (47) 

(=1 l &j i=i 

Here £ v = £oln(A,/£) is the energy of a vortex without a magnetic dot, £o = §l/(16n 2 X) and £ vv 
is the vortex- vortex interaction and £ mv is the vortex-magnetic dot interaction. Substituting d25l) 
into the vortex energy £ vv of (l27t . we get 



£o 
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r ij\ v , T lj 



nj) = -r Ho(£)-Yo(£ 



~2V -ivr (48) 

where ry = |r,- — r ; | , Hq(x) and Yq(x) are the Struve function of the zeroth order and the modi- 
fied Bessel function of the second kind of the zeroth order, respectively [ 59] . For £ mv of d29l . di- 
rect substitution of the vector potential, magnetic field and the phase gradient (see Eqs. (1411441) ) 
into d2"9l gives 

Jo 1 + 2Aq 

In order for N vortices to appear, the necessary condition is that A^ < and A^ < A#-i. 
Using this criteria, we can determine in what configurations and order the vortices appear. To 
this end, we study only vortices with positive vorticity which are situated under the dot. Under 
the assumption that the dot's diameter is larger than \, it is reasonable to think that vortices with 
multiple vorticities, the so-called giant vortices, do not appear, since the vortex energy grows 
as the square of its vorticity (see Eq.d4*71) N ). For large dots R >> X and sufficiently small m, 
the giant vortex is definitely energy unfavorable. However for R < X and mR >> <j>o, it can be 
favorable. This question has not yet been analyzed completely. 

The next step is to minimize d4*7t with respect to the positions of the vortices. We first start 
with one vortex. It turns out that it appears at the center of the dot. Ai is a function of two 
dimensionless variables m§o/e v and R/X. Ai = defines a critical curve separating regions 
with and without vortices, and is depicted in Fig. El Stability occurs (An < 0) for the regions in 
Fig. below the critical curves. The asymptotic behavior of £ mv for large and small values of 
R/X can be found analytically (see Appendix C for details): 

(* » 1), (50) 

- m^ (* « 1). (51) 
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FIG. 3: Magnetic dot on a superconducting film. 



Thus, asymptotically the curve Ai = turns into a horizontal straight line m§o/z v = 1 at large 
RJX and into a logarithmically distorted hyperbola (mfyof&v) (R/X) = 2 at small R/X. 

For further increase of either m§o/E v or R/X, the second vortex becomes energetically favor- 
able. Due to symmetry, the centers of the two vortices are located on a straight line connecting 
the vortices with the center of the dot at equal distances from the center. The curve 2 on Fig. |3] 
corresponds to this second phase transition. The occurrence of 2 vortices can be experimentally 
detected as the violation of circular symmetry of the field. For three vortices, the equilibrium 
configuration is a regular triangle. The further increase of m§o/e v or R/X makes other vortex 
states more energetically favorable. In principle there exists an infinite series of such transi- 
tions. Here, we limit ourselselves to the first seven transitions by considering the next four 
vortex states. In equilibrium, for vortices sit on the corners of a square, whereas five vortices 
form a pentagon. We find that geometrical pattern for six vortices is hexagon. For the case 
of seven vortices, equilibrium configuration is different from the first six states. Namely, one 
vortex is situated at the dot's center, while other six vortices form a hexagonal shape ( see Fig©. 

It is not yet clear what is the role of configurations with several vortices confined within the 
dot region and antivortices outside. Note that these results are valid for infinite systems. When 
the FM dot is placed on top of a finite SC film, we suspect that the geometrical patterns formed 
by vortices will be quite different. To investigate this hypothesis, we recently studied the FM 
dot on a semi-infinite SC film, and found that the vortex is shifted either towards or away from 
the SC film's boundary due to the competition between the attraction of a vortex with the SC 
film's edge via its image vortex and the interaction between the dot and the vortex H26I1 . 
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(a) Two vortices 



(b) Three vortices 



(c) Four vortices 





(d) Five vortices 



(e) Six vortices 



(f) Seven vortices 



FIG. 4: Vortex states 



IV. FERROMAGNETIC-SUPERCONDUCTING BILAYERS 



Earlier Lyuksyutov and Pokrovsky noticed 11201 15311 that in a bilayer consisting of homoge- 
neous superconducting (SC) and ferromagnetic (FM) films with the magnetization normal to 
the plane, SC vortices occur spontaneously in the ground state, even though the magnet does 
not generate a magnetic field in the SC film. In previous work 112 ill , we presented a theory 
of such vortex-generation instability and the resulting vortex structures, and showed that that 
due to this instability, domains with alternating magnetization and vortex directions occur in 
ferromagnetic- superconducting bilayers (FSB). In this chapter, we review the theory of FSB 
and the domain structures. In the next section, we treat these domain structures in the contin- 
uum regime in which the domain size is much larger than the effective penetration depth. This 
approximation does not work when the equilibrium size of the domains is on the order of the 
effective penetration depth. However, it can be recovered by considering the discrete lattice of 
vortices instead. In the second section, we report our preliminary results on the possible equi- 
librium structure in the discrete case and calculate the vortex positions, which depend on the 
magnetization and domain wall energy ||6£ 
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A. The Continuum Regime 

We start by refining previous arguments establishing a topological instability in the FSB 
|20j,|53|]. We assume that the magnetic anisotropy is strong enough to keep the magnetization 



exactly perpendicular to the film (in the z-direction). The homogeneous FM film creates no 
magnetic field outside itself and hence does not alter the state of the SC film. The magnetic 
field generated by a single vortex in the superconducting film, with magnetic flux §q = hc/2e, 
interacts with the magnetization m of the FM film and lowers the total energy by — m§o for 
a proper sign of vorticity. The energy to create a Pearl vortex in an isolated SC film is e v = 
£oln(A,/£) H29I1 . where £o = ^/l6Tl 2 'k, X = X 2 L /d is the effective penetration depth Il57ll . Xl is 
the London penetration depth, and \ is the coherence length. Thus, the total energy of a single 
vortex in the FSB is 

£ e /f = £ v -m^ (52) 

and the FSB becomes unstable with respect to spontaneous formation vortices as soon as £y 
turns negative. Note that close enough to the SC transition temperature T s , ef, is definitely 
negative since the SC electron density n s and, therefore, £ v is zero at T s (Recall that £ V (T) = 

£ v (r = 0)(l — 4j)). For small m value, £ v >0at T = 0, the instability exists in the temperature 

interval T v < T < T s , where Z e v (T v ) — 0, otherwise instability persists until T = 0. 

A newly appearing vortex phase cannot consist only of vortices of one sign. Indeed, any 
system with average vortex density n v would generate a constant magnetic field B z = n v §o 
along the z direction. The energy of this field for a finite film of the linear size Lf grows as 
l}p which quickly exceeds the gain in energy due to creation of vortices, proportional to l£. 
Hence, in order for the vortex array to survive, the film should split in domains with alternating 
magnetization and vortex directions. We show below that if the domain size L is much greater 
than the effective penetration length X, the most favorable arrangement is the stripe domain 
structure. To this end we write the total energy of the bilayer in the form 

U = U sv + U vv + U mv + U mm + Udw, (53) 

where U sv is the effective energy of single vortices; U vv is the vortex-vortex interaction energy; 
U mv is the energy of interaction between the vortices and magnetic field generated by domain 
walls; U mm is the self-interaction energy of the magnetic layer; and Ud w is the linear tension 
energy of magnetic domain walls [60]. We assume that the 2d periodic domain structure con- 
sist of two equivalent sublattices, so that the magnetization m z (r) and density of vortices n(r) 
alternate when crossing from one sublattice to another. The magnetization is assumed to have 
a constant absolute value: m z (r) = ms(r), where s(r) is the periodic step function equal to +1 
at one sublattice and — 1 at the other one. We consider a dilute vortex system where the vortex 
spacing is much larger than X. Then the effective single-vortex energy becomes 

U„ = e e / f J n(r)s(r)d 2 x. (54) 

Note that n(r)s(r) > in all cases. Due to "average neutrality of the periodic stripe system, the 
energy of a single vortex in equation (l54b is different from d52l : £y = £ v — mtyo/2. Note that, 
— \j\s v -md 2 x term in Eq.d2"9t contributes — mtyo/2 in the effective single- vortex energy. In 
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the periodic systems, the contribution of the surface term is zero. The vortex- vortex interaction 
energy is 

U vv = - [n(r)V(r-r')n(r')d 2 xd 2 x', (55) 

where V(r — r') is the pair interaction energy between vortices located at points r and r'. Its 
asymptotic value at large distances | r — r' | ~^> X is H55I1 



V(r-r>) = A2 ^ (56) 

An 1 | r — r I 

This long-range interaction is induced by the magnetic field generated by the Pearl vortices 
and their slowly decaying currents. By Eq.d29l). the energy of the vortex interaction with the 
magnetic field generated by the magnetic film is 12711 



U mv = - T ijL-|vcp(r-r>(r / ) .^{v)d 2 xd 2 x' . (57) 

Here cp(r — r') = arctan ^y is a phase shift created at a point r by a vortex centered at a point 

r' and aS" 1 ' (r) is the value of the vector-potential induced by the FM film upon the SC one. By 
Eq. d28l) . the magnetic self-interaction reads 

Umm = ~ J B { z m \r) S (r)d 2 x. (58) 

Finally, each magnetic domain wall's linear energy is Ud w = EdwLdw, where e^w is the linear 
tension of the magnetic domain wall and Ld w is the total length of the magnetic domain walls. 
Let us analyze the vortex-domain- wall interaction U mv . The magnetic vector-potential A( m ) 
obeys the London-Pearl magneto- static equation (see Eq.(fT51)): 



Vx VxA ( 



4a ( ' n) +47iVx(zm(r)) 8(z). (59) 

K 

We consider L^>X, where the term V x f V x A^ J is negligible and then 

a (m) «-47tm>ixVs(r). (60) 

The phase gradient entering Eq. d57l can be rewritten as: Vcp(r) = £ x Vln |r — r'| . Plugging 
this expression into (1571 . integrating by part and employing relation 

V 2 ln|r — r'| = — 27i8(r — r'), we arrive at 

^0 /".../^../^ ,2 



U mv = -— m(r)n{r)d l x. (61) 

This result implies that the vortex-domain- wall interaction renormalizes the single- vortex to 

e v = £ v -m^ . (62) 

Thus, the term U mv can be removed from the total energy (l53l) if the single- vortex contribution 
U sv is replaced by U v , which differs drom d54b on replacing £y by £ v . In physical terms, 
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it means that the vortex attraction to the magnetic domain walls lowers the threshold for the 
spontaneous appearance of the vortex-domain structure. The next step is the minimization of 
energy with respect to n(r) the vortex density, which appears only in the first three terms of 
the total energy (see Eq.d53"l)). Their sum can be conveniently denoted by U v = U sv + U vv . To 
simplify the minimization, we Fourier-expand the periodic functions: s(r) = Y,G s G elG ' r an d 
n(r) = Y,G n G elG r ■ The energy U v in the Fourier-representation then reads 

U v = Usv + U vv + U mv = £ ( £ v SGn-G + 2^G«G«-G J , (63) 

where Vq = jV(r)e lG ' r d 2 x = <j>q/2tt|G| . Minimization of Eq.d53l over hq leads to 

l v s G 27t£ v |G|s G ., A . 

n G = —rr- = 75 , (64) 

^G <|>o 

U v = -^Z,\G\\s G \ 2 . (65) 

Note that the solution becomes physically meaningless at positive £ v . We now apply these 
general results to analyze the stripe domain structure. In this case the density of vortices n(x) 
depends only on one coordinate x perpendicular to the magnetic domain walls. The vectors 
G are directed along the x-axis. The allowed wave numbers are G = n(2r+ \)/L, where L 
is the domain width and r runs over all integers. The Fourier-transform of the step function 
is sq = w 2 r4-n - The inverse Fourier-transform of Eq. d64l> for the stripe domain case is (see 
Appendix C for details) 

47I£y 1 

«W = — TTT^^E- (66) 

<^Lsmf 

Note the strong singularity of the density near the domain walls. Our approximation is invalid 
at distances of the order of X, and the singularities must be smeared out in a band of the width 
X around the magnetic domain wall. Conversely, the approximation of the zero-width magnetic 
domain wall is invalid within the range of the magnetic domain wall width I. Fortunately, we do 
not need more detailed information on the distribution of vortices in the vicinity of the magnetic 
domain walls. Indeed, by substituting the Fourier-transform of the step function into equation 
(l65t . we find a logarithmically divergent series in the form of £ r 1 / (2r + 1 ) . It must be cut off at 
±>*max with r max ~ L/X. The summation can be performed using the Euler asymptotic formula 
Il65h with the following result (see Appendix C for details): 

4m 2 R / L 



U; T = — I ln- + C + 21n2 1 , (67) 

where m — m — £ v /§o, & is the domain area and C ~ 0.577. Now the problem is to analyze the 
proper cut-off for any lattice. As we have seen in the stripe domain structure, the energy U v 
diverges logarithmically due to the strong singularity of the vortex density near each magnetic 
domain wall (see Eq. d6"6l)). Thus, the logarithmic term is proportional to the total magnetic 
domain wall length (see Eq.doTl)). The singularity of the vortex density contributes the similar 
logarithmic term to the energy for any lattice. However, this logarithmic accuracy is not suf- 
ficient to distinguish the domain structures. In order to solve this problem, we need the next 
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approximation to the energy U v , i.e. a term a, proportional to the length of the magnetic domain 
wall, without the logarithmic factor. Together with this term, the energy for any lattice looks 
like U v ~ ln(La/X). Now, the problem is to find this term accurately. Such a term includes a 
non-local contribution from large distances between X and L and a local contribution from the 
vicinity of the magnetic domain walls. The non-local contribution is accurately accounted for 
by the summation over the integers; whereas, the local contribution requires a cut-off at large r, 
which is not well defined. However, due to its local character it must be the same for all mag- 
netic domain walls. Therefore, it is possible to choose the maximal wave- vector in the direction 
normal to the magnetic domain wall as 2%fk. Such a procedure renormalizes the magnetic 
domain wall's linear tension, in the same for any domain lattice. This remark allows calcula- 
tion of the energy U v for the square and triangular lattices. For the square checkerboard lattice, 
the allowed wave- vectors are G = J [(2r + l)x+ (2s + l)y\. The Fourier-transform of the step 
function is: sq = i(2 ,fw 2 , ]y The maximal values of r and s are identical and equal to Lfk 
where L is the side of a square domain. The summation, similar to the case of stripe structure 
although somewhat more complicated, leads to the following expression (see Appendix C for 
details): 

tf* = -^(ln£+C + 21n2- Y ), (68) 

where the numerical constant y is defined below: 

Y=(2-V5)-5C(3) + 4£ £ S (r,s). (69) 

Here C,(x) is the Riemann zeta-function; ^(3) ~ 1.2020 and, 



z{r A 2(r + ,+ l)- v /(2r+l^ + (2, + l^ 

H,) ~ (2r+ l) 2 (2s + 1) 2 ' { > 

Direct numerical calculation gives y ~ 0.9 > In 2. 

The reciprocal lattice vectors for the regular triangular domain lattice are 



G=f 



^ A "-73>V + *7^ 



The analysis is remarkably simplified in the "triangular coor- 
dinate frame": u = x — y/V3; v = 2y/y3. The step function inside one elementary cell is 
s(u, v) = + 1 for u + v < L and s(w, v) = — 1 for u + v > L, where L is the side of the elementary 
triangle. The Fourier-transform of the step function sq is not zero at either r ^ 0, s = 0, or 
r = 0, s 7^ 0, or r = s ^ 0. For all these cases \sq\ 2 = l/(n 2 q 2 ), where q is either r or s, depend- 
ing on which of these numbers differs from zero. For this case, the summation in equation d65l) 
gives 

t/r = -^^(lnr max +C). (71) 

However, the value r max is different from the stripe and square cases since the coordinates are 

skewed, here it equals ^^. 

Our next step is to show that the magnetization self-interaction can be included into the 
renormalized magnetic domain wall linear tension. For the isolated FM stripe domain struc- 
ture, the magnetization self- interaction energy is equal to U mm = —m 2 L^ vv In j, where I is the 
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magnetic domain wall width 116 U 1621 16311 . The superconducting screening enhances the mag- 
netic field very near the magnetic domain walls and reduces it way from the magnetic domain 

walls. In the stripe geometry, from Eq. d60l) . b" = daS m > /dx = —A%mk(d 2 s/dx 2 ) implying 
that the screened magnetic field is confined to an interval ~ X near magnetic domain walls. 
Thus, its contribution to the energy does not contain a large logarithm. By (l58l) . U mm gives 
—Ld w m 2 , which can be incorporated into the renormalized value of the magnetic domain wall 
linear tension. Note that this contribution is negative. We assume that it is less than the initial 
positive linear tension £^ w . We do not consider here the interesting but less likely possibility of 
a negative renormalized linear tension, which probably results in domain wall branching. 

Now we are in position to minimize the total energy U over the domain width L and compare 
the equilibrium energy. The equilibrium domain width and the equilibrium energy for the stripe 
structure are 

(str) ^ ( Edn 



L^ = ^exp(^-C+l), (72) 



e> = -^exp(-^ + C-l). (73) 

Calculating the corresponding values for the square and triangular lattice, we obtain: L e s q — 

Leg exp(y); Ueq = 2Ueq exp(— y); and Ueq = (3/4) {/j^. Comparing these energies to 

the energy of the stripe structure, for which Ueq < 0, we conclude that the stripe structure wins. 
The domains become infinitely wide at T = T s and at T = T v . The expression in the exponent 
( 1731 at T = T s is four times less than the corresponding expression for domains in an isolated 
magnetic film 15611 . Therefore, stripe domains in the bilayer can be energetically favorable 
even if the isolated magnetic film remains in a monodomain state. If stripe domains in the 
magnetic film exist above the SC transition, then they shrink dramatically below the transition. 
The physical reason behind this effect is as follows. There are two contributions to the energy 
that determines the domain width; magnetic energy and the domain wall energy. While the 
latter contribution prefers the larger domain width, the former contribution causes the smaller 
width. In the presence of superconductor, the magnetic field is screened due to the Meissner 
effect, which makes the magnetic energy contribution dominant. As a result, the domain width 
is smaller than that in the case of ferromagnet only. Bulaevsky and Chudnovsky [36] found that 

1/3 

the domain width in a thick magnetic layer above a bulk superconductor is proportional to dm 

1/2 

instead of d m , a well-known result for an isolated magnetic layer. Here d m is the thickness of 
the magnetic layer; "thick" means that d m »l. Our problem is fundamentally different on two 
counts: first, we consider a thin FM film d m <C I above a thin SC film and, second, the main 
effect is due to the interaction of vortices with the magnetization rather than from the screening 
of the magnetic field as in IB6I1 . The vortex-magnetization interaction effect is much stronger, 
leading to a totally different dependence. 



B. The Discrete Regime 

If Erfw < 4m 2 , the continuum approximation becomes invalid, since L eq becomes on the order 
of or less than X (see Eq.(|72l). Instead a lattice of discrete vortices must be considered. In this 
section, we present a method which works in both continuum and the discrete regimes. We study 
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FIG. 5: The vortex lattice. 



the lattices of discrete vortices only in the stripe phase, since the system favors that phase in 
equilibrium. In the continuum approximation, it is found that the vortex density increases at the 
closer distances to the magnetic domain walls. Based on this fact and the symmetry of the stripe 
domain structure, it is reasonable to consider that the vortices and antivortices form periodic 
structures on straight chains along the y direction. Even though it is not clear how many chains 
are associated with each domain, we can still make progress toward understanding discrete 
vortex lattices. In this section, we report our prelimanary results on the problems; i) how the 
vortices and the antivortices are positioned on the chains; ii) how the equilibrium domain size 
changes, depending on the magnetization and the magnetic domain wall energy in the presence 
of the vortices. In order to solve these problems, we propose a configuration of the vortex and 
the antivortex chains in which two chains per stripe is considered. According to our results in 
the continuum regime, the vortex density increases near the domain walls. Therefore, we can 
have at least two chains per stripe in a possible configuration. Chains are situated at distance a 
from the magnetic domain walls. Another problem is how to place vortices and antivortices on 
the chains. If vortices are next to each other on the either side of the magnetic domain wall, the 
magnetic fields they produce cancel out each other. As a result, the gain in energy is diminished. 
However, the system has the largest gain when the vortices and antivortices are shifted by a half 
period b/2 (see Fig. |5]). 

Our next step is to write the energies of the proposed structure. To this end, we use 
(I36I37I38I ). In those equations, the vortex configurations differ by their form-factors. We can 
obtain them from Fq = £ r .M,-e' G ' r ', where the G's are the reciprocal vectors of the periodic 
structures, the r ; are the positions of the vortex centers, and n,- are the charge of the vortex. In 
our proposed model, G= ((2r+l)f,2sf),w,- = ±l andF G = e iG * a - {-l) s e~ iGxa . Inourcalcu- 
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7i(2r+l) 



lations, the divergent part of the series must be extracted carefully. We show below the detailed 
analysis of series equations for each candidate. We start with the self interaction energy of the 
magnetic layer U mm , since it is the same for each configuration. For the periodic structures, it 
is given by ( BHl . Direct substitution of the Fourier coefficient of the stripe phase m z Q 
into Eq.d^St gives the self-interaction of the magnetic layer per unit cell as 

8m2 V ! 

L toik + ^+l 

This series is logarithmically divergent. However, it can be pulled out easily by adding and 
substracting l/(2r+ 1) in the series above. Thus, we get two terms; one convergent, the other 
divergent. Summing over r on the divergent part up to the cutoff r max = L/l, where / is the 
magnetic domain wall width, we obtain the following: 

«- = -£( ta 7-'K + i3C>)' (?5) 

where \|/°) (x) is the polygamma function of zeroth order [ 59] . In our numerical calculations, we 
write the logarithmic term in (1751 as ln(A,/Z) + ln(L/X) and then incorporate the — 4m 2 ln(X//) 
term in the renormalized eJ" Another energy term with a divergent series is the vortex energy, 
in general given by (BBT) . The logarithmic divergence in this term stems from the vortex self- 
energies. We first split (I3BT ) into two parts as follows: 



_rt£0_y 

2L 2 b 2 % 



\Fg\ 2 \F G \ 



G 2 G 2 (l+2XG) 



(76) 



Note that the area of the elementary cell is 2Lb. The first term of the series above contributes 
to the self-energies of the vortices; whereas, the second term is the vortex- vortex energy and 
will be left in the series form. The series in the first term can be transformed to the form 
of L^L-ooLIL-oo l/((2f+ l) 2 x 2 + s 2 ) where x is constant, and depends on the form-factor. A 
detailed analysis of such series is given in Appendix B. 

The next step is to find the vortex energy and the interaction energy of the magnetization and 
vortices for each configuration. In the calculation of u mv , we take the Fourier coefficient of the 
magnetization to be nr+u ^i^y)- ^ ne fact that the stripe is infinite along the y direction results 
in the additional term 2%8(G y ). However, it does not play any role in the calculation of u mm . 
For numerical analysis, these energies must be expressed in terms of dimensionless parameters. 
To this end, we define dimensionless variables X = X/L , b = b/L and e^w = e rfw^/ e °- ^ ne tota ^ 
energy U is measured in units of £q/X 2 . In addition, we introduce the dimensionless magnetic 
energy as U mm = u mm /(&o/X ). In the fourth configuration, the square of the form- factor is: 
I^FqI 2 = 2 — 2( — l) i? cos((2r + l)nd). Even and odd values of s give different contributions. 
Then, we can calculate the vortex energy for even s and odd s separately. In terms of these 
parameters, we find 

U = K fln(4) -2f v (X,a) - ^-f vv (X,a,b) - ^^f mv (X)) +u mm + £ dw X 1 (77) 

2b V XZ, bn £o / 
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where 



_ - coth((2r+l)f) 
Jv — ]_, 

r=0 



2r+l 



■sin 2 ((2r+l)7Ta) 



+ 1 

r=0 



tanh((2r+l)f)-l . 

— sin ((2r+l)7ta) 



2r+l 



/v- 



sin 2 ((2r+l)7Ta) 



r^-((2r + l)2+^)(l+27dLJ(2r+l)2 + ^) 



+ I 

r,s=— < 
/mi' — / , 



cos 2 ((2r+l)rca) 



4(2.-+ 1) 2 



((2r+l) 2 +^^)(l+2^ v /(2r+l) 2 + 
sin((2r+l)7la) 



4(2.y+l) 2 



^ (2r+l)(l+27r?i(2r+l)) 



(78) 



In the numerical minimization of Eq. dTTl . we take ln(A,/^) = 5. Changing m^>o/ e O at fixed 
e^ w , which initially is fixed at 0.5, we calculate the minimal energy of the proposed config- 
uration. We first investigate when this configurations becomes energetically favorable in the 
system. To this end, we check where the equilibrium energies of the configuration first become 
negative. 

TABLE I: The numerical results for the fourth configuration at £^ vv , = 0.5. The column on the left is 
input. 



H>o/e v 


L/X 


a/k 


b/X 


2.00 


2.50 


0.56 


0.18 


3.00 


1.67 


0.38 


0.12 



In numerical calculations, we also found that the vortex lattice is stable for e^w > 4m 2 . At 
this point, the domain size is noticably larger than the effective penetration depth A,, so the 
continuum approximation is valid. Therefore, we expect that the domain nucleation starts in the 
continuum regime. This problem is left for the future research. As seen in Table |IJ at constant 
td w , with increasing m§o/£ v , the equilibrium size of the domain decreases. In addition, the 
vortices on the chain get closer to each other. These results agree with those obtained in the 
continuum approximation. As £rf w /4m 2 increases, we expect that new vortex chains develop 
within the domains. We leave this problem to future research. 

V. CONCLUSIONS 

We reviewed theory of the heterogeneous magnetic superconducting systems (HMSS) based 
on London-Maxwell equations and the application of the theory on two realizations: ferro- 
magnetic (FM) dots and their square array on a superconducting (SC) film, and ferromagnetic- 
superconducting bilayers (FSB). In the first chapter, we presented a general formalism for the 
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interaction between magnetic textures and superconductors in the Londons approximation. The 
problem is formulated as a variational principle. The variational functional (energy) is an inte- 
gral over regions occupied either by magnet or by superconductor. It allows us to find directly 
the positions of vortices and magnetization. Afterwards, the formalism is extended to the case 
of periodic structures and finite systems. 

As applications of the formalism, we have shown that vortices in superconducting films can 
be generated by magnetic dots magnetized normal to the film. We have found phase transition 
curves separating the state without vortices from the state with one vortex and the latter from 
the state with two vortices. Up to 7 vortices, the vortex-configurations on the ground state of 
a SC film with the FM dot on it are determined. For one vortex under a dot we have shown 
that the perpendicular component of magnetic field changes sign at some distance from the dot. 
This fact can be used for diagnostics of the vortex generation. 

However, we treat only vortices under the FM dot. In a more realistic picture, the antivor- 
tices outside the dot become important and most likely affect the configurations of the vortices 
confined within the dot's region. This problem still remains open. 

In the fourth chapter, we studied ferromagnetic- superconducting bilayers (FSB). We pre- 
dicted that in a finite temperature interval below the SC transition the FSB is unstable with 
respect to SC vortex formation. The slow decay (<x \jr) of the long-range interactions between 
Pearl vortices makes the structure that consists of alternating domains with opposite magneti- 
zation and vorticity energetically favorable. The distribution of vortices inside each domain is 
highly inhomogeneous, with density increasing near the magnetic domain walls. As long as 
the domain width is larger than the effective penetration depth, the energy of the stripe domain 
structure is minimized. These new topological structures can be observed directly. A strong 
anisotropy in current transport would provide indirect evidence of the stripe texture: the bilayer 
may be superconducting for current parallel to the domains and resistive for current perpendic- 
ular to the domains. 

If Erfw < 4m 2 , the continuum approximation becomes invalid. Instead, we considered the dis- 
crete lattice of vortices. We analyzed the vortex configurations in two vortex chains. Depending 
on the magnetization and the magnetic domain wall energy, the positions of the vortices and the 
equilbrium domain size are calculated. 

It is possible that the long domain nucleation time can interfere with the observation of 
described textures. We expect, however, that the vortices that appear first will reduce the barriers 
for domain walls and, subsequently, expedite domain nucleation. Quantitative study of this 
dynamic process is still in progress. 

Our purpose was to consider quantitatively a new class of phenomena provided by the inter- 
action between superconductivity and ferromagnetism in heterogeneous systems. This review 
article focuses on unusual equilibrium structures of vortices and magnetization occuring in such 
systems. As the simplest specific examples we considered single magnetic dots and also SC- 
FM bilayers. We have shown that, because of the numerous parameters for the system, such as 
temperature, magnetization, thickness etc., these systems display a complex phase diagram. We 
did not exhaust all possible states. Moreover, dynamic effects and transport properties are be- 
yond the scope of this article. However, we believe that our primary results are of experimental 
interest and have technological promise. 

We believe that the most experimentally interesting and challenging predictions are: (1), the 
existence of domain structures with alternating magnetization and vortex polarity in FM-SC 
bilayers and their shrinking below the SC transition temperature; (2), the formation of several 
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vortices under the magnetic dots placed on top of a SC film; (3), the symmetry violation in 
the periodic array of the magnetic dots on SC films. Although, antivortices outside the single 
magnetic dot are not analyzed in this article, they are expected to occur in real systems. All 
these effects can be observed directly by scanning tunnelling microscopy, scanning Hall probe 
microscopy and micro SQUID measurements. 

In conclusion, our results not only confirm some old results found by means of different 
methods, but also present a class of new physical effects in HMSS. In this sense, they manifest 
a new direction and motivation for the possible experiments in the future. 

APPENDIX A: THE PEARL VORTEX 

The vortices in thin superconducting(SC) films are first studied by J. Pearl 12911 . Here, we 
give the detailed calculations of vector potential and magnetic field for Pearl vortex located at 
ro on the SC film. We start with the London-Pearl equation: 

V 2 A(r,z) = I a (r)5(z) -^V(p(r-r )5(z), (Al) 

where a(r) = A(r,z = 0). It is easy to find the vector potential due to the vortex by employing 
the Fourier transformation to Eq. dAll) . In doing so, we use 

A(,,r) = /A**«^. (A2) 

In the Fourier representation, the London-Pearl equation reads 

1 4> (Vcp)q rA ~ 



where (Vcp) q = 2n l -^^e ,q ' r ° 115 711 . Employing integral over k z to Eq. (IA3l) and using a q = 

I-<=° ~2k ^ k ' we °btain the Fourier-transform of the vortex-induced vector potential at the SC 
film as 



%)(<? x z) g -'q - r o 

q{\ + 2Xq) 
Substituting the above equation in Eq. (IA3l) . the 3d vortex-induced vector potential is found as 



a q = zk , k In ■ ( A4 ) 



2iUgxz)e-^ 
Ak= k2(l+2^) • (A5) 



The direct substitution of Eq. (IA5l) back in Eq. (IA2l) leads to 



2iU^z)e-'^-^^ d 2 q dk z 
A(r ' ZJ " ' k2(l+2M 72^' (A6) 



First, we perform integral over k z , and find 



4o(^l)^Mr*l d 2 q 
A(r ' - ] I q(l+Zkq) (2J0 2- C } 
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Let 



q = gcos(9 + (|))w + gsin(0 + (|))v (A8) 

r — r = |r — r |cos(|)a + |r — r |sin(|)v, (A9) 

where u and v are unit vectors in a plane perpendicular to the z-direction, 9 is the angle between 
q and r — ro, and § is the angle between r — ro and u. By (IA8I) . q x z in Eq. (IA7l) reads 

<?x| = -cos(0 + 4))v + sin(0 + 4))M. (A10) 

Substituting dATOb in Eq.dA7l>. and using u x | = -v, v x z = u, f^ e ixcosQ cosQdQ = 2%iJi(x) 
and J^ n e ixcosB sinGJG = 0, we find 

= ^-x(r-ro) f-y.tolr-rol).-*! 

271 I r — ro I ./o l+2Ag 



Tx(r-r ) 
components of the Pearl vortex are found as follows: 



Note that in the above equation, cos(|)v — sintyii = i _ r | . From B = V x A, the magnetic field 



flr(r,z) = ^stgn(z)y o fT ^ *Z, (A12) 

5z(r ' z) = 2777o TT2^ dq - (A13) 

APPENDIX B: INTEGRALS OF BESSEL FUNCTIONS 

In this appendix, the asymptotic values and the exact results of the Bessel integrals used in this 

article are introduced. First, we present the integrals containing one Bessel function. These 
integrals are in the form of 

' Ukr)kn -dk, (Bl) 



/o 1+2J& 
where m = 0, 1 and n = 0, 1 . For m = 0, 1 and n = 0, the exact results can be obtained as follows 



[64]: 



o 1+2U 4A L y 2X J uv 2r J ' 

o 1+2J& 4A L U 2A/ V 2A ;J r 

where i//(jc) and F/(x) are the Struve function and the second kind of the Bessel function of the 

/th order [59]. The respective asymptotic values of Eqs. (IB2[|B3l) are found as follows: When 



r « A, which corresponds to kX >> 1, the integral in dB2l> becomes [65] 



„ ^*=aN7>-<* 
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For r » X, which is equivalent to kX« 1, the integral in dB2l> becomes W65I1 



J (kr)dk=-. (B5) 

o r 



Using the same techniques, the asymptotic values of the integral in dB3l) are 



Hkr)_ dkK rm» = » r<<l ( B6) 



o 1 + 2kX Jo 2kX 2X 

Jl{kf > dk re / Ji(kr)dk = - r»X. (B7) 



/o 1 + 2&X 7o a- 

In order to find the asymptotic values for m = 0, 1 and n = 1, we use 

/ J m {kr)k n dk = Yim a ^ Q — J m (kr)e- ak dk. (B8) 

Jo oor Jo 

With similar techniques and (IB8b . the asymptotic values are given as follows: 

7 °(* r Wrm« = j- r «x. <B9) 



o 1 + 2fcX Jo 2X 2Xr 



For r » X, the integral becomes [q Jo(kr)kdk, which equals zero [65]. In order to find nonzero 
result, we do the following approximation: r » X is equivalent to kX « 1. Therefore, the 
fraction in (IB II) can be rewritten as 

1 \-2kX+--- (BIO) 



1+2U 
By (IBlQb . we obtain 



/ J °( kr { kdk re -2A, / J (kr)k 2 dk=^ T r»X, (Bll) 

Jo 1 + 2£A Jo r J 



./o 
and, 



*(*>,«»„ rfa = -J. r «x ( Bi2) 



1 + IkX Jo IX 2Xr 

/ „ u ' kdk re / J l (kr)kdk=^ r»X. (B13) 

Jo 1 + 2&A Jo r z 

We give only the asymptotic values of the integrals with two Bessel functions. We start with 
the following integral: 

Jo l+2kX dL (BU) 

For the above integral, we first analyze the case, in which R « X. In this case, we can replace 
7i (kR) by kR/2 in (IB 141) . In doing so, we get 

f" Ji{kR)Jo{kr) R r Mkr) 

Jo l+2kX dk ~2j TT2kX kdL (B15) 
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Using the asymptotic values for the integral on the left in (IB9I ) and (IBllL we find 



• Ji(kR)J (kr) R 

1 + 2kX 4Xr 

1X.R 
= s- R«X«r. (B17) 

r i 



For R > X and r > A, we can neglect 2kX in Eq. (IB14l) . In doing so, we obtain 

' Ji(kR)J (kr) 



f™P&« f Ji{kR)J {kr)dk. (B18) 

Jo 1 + 2fcA Jo 



The above integral equals [65] 



/ Ji(kR)J (kr)dk = R < r (B19) 

Jo 

= ^ * = ' (B20) 

r<R. (B21) 



7? 
The other integral of interest containing two Bessel functions is 

Jx{kR)J Q {kr)k 2 



/o l+2kX 

For R <X, using 7i (fc/?) ps fci?/2, the integral in Eq. (IB22l) becomes 



dk. (B22) 



R rwr)e_ dk 



2 Jo l+2kX 

By (IB 81) and (IBIOL the asymptotic values of the above integral can be calculated as follows: for 
r < < X, (IB22I) can be rewritten as 

— / J (kr)k 2 dk = — 1 . (B24) 

4A Jo 4Ar J 

For R << r, (IB23I) becomes 

-i?A / J (kr)k 4 dk= =-. (B25) 

Jo 7° 



Now, we can write the asymptotic values of dB22l) as 

r Ji(kR)J (kr)k 2 R , /T ,_ S 

r J -Wy^-^. *«X«r. (B27) 

Jo 1+ 2a:A r 3 
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APPENDIX C: CALCULATION OF SERIES 

In this appendix, the detailed analysis of series is given. First, the series in the energy cal- 
culations of the periodic systems are analyzed; second, the detailed calculation of the vortex 
density is shown. The series we encounter in the energy calculations fall into two categories. 
In the first category, we sum over one variable. The series in this category are in the form of 
L^=i l/ r - Employing the Euler-Maclaurin summation formula [66], the summation is found 
with logarithmic accuracy as 

£-~lnw + C. (CI) 

where C ~ 0.577 is the Euler-Mascheroni constant. If the summation is performed over only 
odd integers, we can still transform our series to (IC1I) . Namely, 

r max 1 2r mflv +l i i r max /2 , 

« ln(2r mflX + l)+C-ln(^)-^, (C3) 

« ^(lnr max + C + 2\n2). (C4) 

The second category is the double series. In this aspect, we first show the calculation of square 
domain energy in the continuum approximation. The corresponding energy contains the series 

- y/(2r+l)2 + (2j+l) 2 (C r, 

r to (2r + 1)2(2*+ 1)2 • CU>J 

Our goal is to calculate the logarithmic contribution due to the self-vortex energies and the 
constant as a next approximation. The sum in dC5l ) diverges logarithmically in two regions: 
r » s and s >> r. Keeping this in mind, (IC5I ) can be rewritten as 

^ max(2r+ 1,25+1) 
" rj k. (2r+ 1)2(25+ 1)2 



f- V(2r+ 1)2 + (2s+ l) 2 -max(2r+ 1,2^ + 1) 
! L (2r+ 1)2(25 + 1)2 • (C0J 



r,s=— oo 



The first term in the above series contributes the logarithmic term S/ og , and here it is 



s, " s = 2 L^to(^w (C7) 

~ — (lnr mflX + C + 21n2). (C8) 
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We used the result in Eq. (IC4l) . and £^L jtttv. = T ^^ - Th e otner term in (IC6I) contributes 
the constant S C0MiS . For three regions r = s, r > s and s > r, the second series in (IC6I ) is rewritten 
as follows: 

S_-(V2 2)^ (2r+l)3+ 2^^ (2r+ 1)2(2,+ 1)2 " (C9) 

Numerical calculation gives S CO ns = — 1.19. 

The other double series of interest here are in the form of 



f— oo t y=oo 



'W = I E To' <C10) 



7"= — 005=: — 00 " 



where x is an arbitrary constant. Although (IC10I) is logarithmically divergent, the sum over one 
of the variables can be done easily. To this end we perform the sum over s first. In doing so, 



Eq. (IC10I ) becomes (2n/x)J^LiCoih(Tixr)/r JI67I1 . This series is logarithmically divergent. In 



order to get the logarithmic term , we add and subtract \/r. Using the result in (IC1I) . finally we 
get 

2n r ^ coth(rocr) - 1 
I(x) « — E Hhlnr^ + C]. (Cll) 

* r= 1 ' 

Employing the same techniques, we give the results of the different versions of Eq. (IC10I) 
below: 



r f° y 1 2ji " coth((2r+l)jtt)-l ln W C 

i^A* 2 (2r+l) 2 + * 2 ~ *^ 2r+l + 2 + 2 J ' ^ 

Y S y 1 7t - tanh((2r+l)f)-l \nr max C 

i^i^ 2 (2r+l)2 + (2. + l)2~*12- 2r+l ~ r 2 + 2 J " ^ ^ 

In (ICT2l and (ICT3T) . we use £~ =0 l/(v 2 + (2, + l) 2 ) = 7Uanh(7ly/2)/(4y). In the presence of 
sin 2 ((2r+l)y) orcos 2 ((2r+l)v), the series can be calculated in a similar way, using sin 2 ((2r+ 
\) y ) = (l-cos(2(2r+l)y))/2orcos 2 ((2r+l)v) = (1 +cos(2(2r+ l)y))/2. For example, 

r y s y sin 2 ((2r+l)y) _ 2tc " sin 2 ((2r+ l)y)(coth((2r+ 1)jw) - 1) 
2- 2- (jifo r -L-\\2 _i_ a\ ~ v 12- 



/•= — 00 ,y= — co 



(x 2 (2r+l) 2 + 5 2 ) jc \% 2r+l 



+ lnr max ln|cot(y/2)| | C 



r y™ s y cos 2 ((2r+l)y) _ 2% " sin 2 ((2r+ l)y)(coth((2r + 1)jw) - F 



/' = CO £= CO 



(x 2 (2r+l) 2 + s 2 ) x L r to 2r+l 



lnr max In I cot(y/2) I C n 

_j^ + J — ^v_^i+ . (C15) 

4 4 4 J 
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Wb use 

f, cos((2r+l)9) ln|cot(9/2)| 

£ 2r+l = 2 • (C16) 

The Fourier transform of the vortex density hq for the stripe domain structure is found by 
substituting the corresponding reciprocal lattice vectors G = n(2r+ \)/L and the Fourier trans- 
form of the step function sq = ^r+n m Eq-(R>4"1>- In doing so, we obtain 

nG = -TS- I ^n(2r+l). (C17) 

Employing the inverse Fourier transform to Eq. dC17l) gives the following series: 

oo oo 

F(x) = £ sign(2r+l)e'( 2r+1 )T = 2/£sin((2r+l)^). (C18) 

In order to calculate the above series, we integrate both sides over x. In doing so, by dC16l ) we 
find 

f , s ^ cos((2r+l)^) L , ,«c„ 

/ F *)<& = -2i- £ " , /^ = -/-In cot - • (C19) 

J 71 £? 2r + 1 7T 2L 

The derivative of the above equation with respect to x gives 

m = --^y- (C20) 



By (IC17I) and (IC20I) . the vortex density m(jc) becomes 

47t£ v 

WoMf 



n{x) = -^-^=r. (C21) 



30 



IV. REFERENCES 



[1 

[2 



[3 

[4 

[5 
[6 
[7 
[8 
[9 
[10 

[11 
[12 
[13 

[14 

[15 
[16 
[17 
[18 
[19 

[2o; 

[21 

[22 
[23 
[24; 
[25 
[26 
[27 
[28 
[29 

[30; 

[31 
[32 
[33 



I.F. Lyuksyutov and V.L. Pokrovsky, Phys. Rev. Lett. 81, (1998) 2344. 

I.F. Lyuksyutov and V.L. Pokrovsky, Superconducting Superlattices II: Native and Artificial, Vol. 

3480 (Eds. Ivan Bozovic and Davor Pavuna), SPIE-International Society for Optical Engineering, 

Bellingham, WA, (1998), pp. 230. 

L.N. Bulaevskii, A.I. Buzdin, M.L. Kulic and S.V Panyukov, Adv. Phys. 34, (1985) 175. 

0. Fisher, Magnetic Superconductors, Vol. 5 (Eds. K.H.J. Buschow and E.P Wohlfarth), North- 
Holland, Amsterdam, (1990), pp. 465. 

M. Belie and M.V Jaric, Phys. Rev. Lett. 42, (1979) 1015. 

H.S. Greenside, E.I. Blount and CM. Varma, Phys. Rev. Lett. 46,(1981) 49. 

V.L. Ginzburg, Zh. E'ksp. Teor. Fiz. 31, (1957) 202 [Sov. Phys. JETP 4,(1957) 153]. 

PW. Anderson and H. Suhl, Phys. Rev. 116, (1959) 898. 

P. Fulde and R.A. Ferrell, Phys. Rev. 135, (1964) 550. 

A.I. Larkin and Yu.N. Ovchinikov, Zh. E'ksp. Teor. Fiz. 47, (1964) 1136 [Sov. Phys. JETP 20, 

(1965) 762]. 

H. Suhl, J. Less-Common Metals 62, (1978) 225. 

B.T. Mathias, H. Suhl and E. Corenzwit, Phys. Rev. Lett. 1, (1958) 92. 

C. Bernhard, J.L. Tallon, Ch. Niedermayer, Th. Blasius, A. Golnik, E. Brucher, R.K. Kremer, D.R. 

Noakes, C.E. Stronach and E.J. Ansaldo, Phys. Rev. B 59, (1999) 14099. 

L. Bauernfeind, W. Widder and H. Braun, Fourth Euro Ceramics, Vol. 6 Faenza Editrice S.p.A., 

Italy, (1995). 

1. Felner, U. Asaf, Y. Levi and O. Millo, Phys. Rev. B55, (1997) R3374. 
S. Erdin, Physica C 391,(2003) 140. 

I.K. Marmorkos, A. Matulis and FM. Peeters, Phys. Rev. B 53,(1996) 2677. 

M. A. Kayali, Phys. Lett. A, 298, (2002) 432. 

M. A. Kayali, Phys. Rev. B, 69, (2004) 012505. 

I.F. Lyuksyutov and V.L. Pokrovsky, cond-mat/9903312 (unpublished). 

S.Erdin, I.F. Lyuksyutov, V.L. Pokrovsky and V.M. Vinokur, Phys. Rev. Lett. 88, (2002) 017001. 

M.A. Kayali and V.L. Pokrovsky, Phys. Rev. B 69, (2004) 132501. 

V.L. Pokrovsky and H. Wei, Phys. Rev. B 69, (2004) 104530. 

I.F. Lyuksyutov and D.G. Naugle, Modern Phys. Lett. B 13, (1999) 491. 

I.F Lyuksyutov and D.G. Naugle, Int. J. Mod. Phys. B 17, (2003) 3713. 

S. Erdin, in press Phys. Rev. B 69 (2004). 

S.Erdin, A.M. Kayali, I.F. Lyuksyutov and V.L. Pokrovsky, Phys.Rev.B 66, (2002) 014414. 

E.H. Brandt, Rep. Prog. Phys. 58, (1995) 1465. 

J. Pearl, Appl. Phys. Lett. 5, (1964) 65. 

M.V Milosovic and FM. Peeters, Phys. Rev. B 69, (2004) 104522. 

J.E. Santos, E. Frey and F Schwabl, Phys. Rev. B 63 (2001) 4439. 

R. Sasik and T. Hwa, cond-mat/0003462 (unpublished). 

G. Carneiro, Physica C 404, (2004) 78. 



31 



[34] D.J. Priour Jr. and H.A. Fertig, Physica C 404, (2004) 293. 

[35] L.N. Bulaevskii, E.M. Chudnovsky and M.P Maley, Appl. Phys. Lett. 76, (2000) 2594. 

[36] L.N. Bulaevskii and E.M. Chudnovsky, Phys. Rev. B 63 (2001) 2502. 

[37] L.E. Helseth, Phys. Lett. A, Vol.315, (2003) 399. 

[38] J.I. Martin, M. Velez, J. Nogues and I.K. Schuller, Phys.Rev. Lett. 79, (1997) 1929. 

[39] D.J. Morgan and J.B. Ketterson Phys. Rev. Lett. 80, (1998) 3614. 

[40] Y. Otani, B. Pannetier, J.P Nozieres and D. Givord, J. Magn. Mag. Mat. 126, (1993) 622. 

[41] O. Geoffrey, D. Givord, Y. Otani, B. Pannetier and F. Ossart, J. Magn. Mag. Met. 121, (1993) 223. 

[42] Y. Nozaki, Y Otani, K. Runge, H. Miyajima, B. PannetierJ.P Nozieres and G. Fillion, J. Appl. 

Phys. 79, (1996) 8571. 
[43] M.J. Van Bael, L. Van Look, K. Temst, M. Lange, J. Bekaert, U. May, G. Guntherodt, V.V. 

Moshchalkov and Y Bruynseraede, Physica C 332, (2000) 12. 
[44] A. Terentiev, D.B. Watkins, L.E. De Long, D.J. Morgan andJ.B. Ketterson, Physica C 332, (2000) 

5. 
[45] M.J. Van Bael, K. Temst, V.V. Moshchalkov and Y Bruynseraede, Phys. Rev. B 59, (1999) 14674. 
[46] PMartinoli, N. Nsabimana, G.A. Racine, H. Beck and R. Clem, Helv. Phys. Acta. 56, (1983) 765. 
[47] V.V. Metlushko, M. Baert, R. Jonckheere, V.V. Moshchalkov and Y Bruynseraede, Solid State 

Comm. 91, (1994) 331. 
[48] M. Baert, V.V. Metlushko, R. Jonckheere, V.V. Moshchalkov and Y Bruynseraede Phys. Rev. Lett. 

74, (1995) 3269. 
[49] V.V. Moschalkov, M. Baert, V.V. Metlushko, E. Rosseel, M.J. VanBael, K. Temst, R. Jonckheere 

and Y Bruynseraede Phys. Rev. B 54, (1996) 7385. 
[50] V.V. Metlushko, L.E. DeLong, M. Baert, E. Rosseel, M.J. VanBael, K. Temst, R. Jonckheere and 

Y Bruynseraede, Europhys. Lett. 41, (1998) 333. 
[51] I.F. Lyuksyutov and D.G. Naugle, Physica C 341, (2000) 1267. 
[52] I.F. Lyuksyutov, D.G. Naugle and V.L. Pokrovsky, Superconducting and Related Oxides: Physics 

and Nanoengineering IV, Vol. 4058 (Eds. Davor Pavuna and Ivan Bozovic), SPIE-International 

Society for Optical Engineering, Bellingham, WA, (2000), pp. 376. 
[53] I.F. Lyuksyutov and V.L. Pokrovsky, Modern. Phys. Lett. 14, (2000) 409. 
[54] L.E. Helseth, PE. Goa, H. Hauglin, M. Baziljevich and T Johansen, Phys. Rev. B 65, (2002) 

132514. 
[55] PG. de Gennes, Superconductivity of Metals and Alloys, Addison-Wesley, (1989). 
[56] L.D. Landau and E.M. Lifshitz, Electrodynamics of Continuous Media, Butterworth-Heinemann, 

(1995), 2nd ed. 
[57] A.A. Abrikosov, Introduction to the Theory of Metals North Holland, (1986). 
[58] V.G. Kogan, Phys. Rev. B 49, (1994) 15874. 
[59] M. Abramowitz and I. A. Stegun, Handbook of Mathematical Functions, Dover Publications, 

(1970). 
[60] A.P Malozemoff and J.C. Slonczewski, Magnetic Domain Walls in Bubble Materials, Acedemic 

Press, (1979). 
[61] Y. Yafet and E.M.Gyorgy, Phys. Rev. B 38, (1988) 9145. 
[62] A. Kashuba and V.L. Pokrovsky, Phys. Rev. B 48, (1993) 10335. 
[63] A. Abanov, V Kalatsky, V.L. Pokrovsky, W.M. Saslow, Phys. Rev. B 51, (1995) 1023. 
[64] YL. Luke, Integrals ofBessel Functions, McGraw-Hill, (1962). 



32 



[65] I.S. Gradstein and I.M. Ryzhik, Tables of Integrals, Sums, Series and Products, Academic Press, 

(1994), 5th ed. 
[66] G.B. Arfken and HJ. Weber, Mathematical Methods for Physicists, Academic Press, (2000), 5th 

ed. 
[67] E.R. Hansen, A Table of Series and Products, Prentice-Hall, (1975). 



33 



